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Abstract 

We report a dynamical formation of a (Mott) gap from holographic fermions. By coupling a 
fermion field with dipole action to the charged dilaton black branes with a Lifshitz like IR geometry 
and AdSi boundary, we find that when the dipole interaction is large enough, spectral weight is 
transferred between bands, and beyond a critical dipole interaction, a gap emerges in the fermion 
density of states. The value of the gap becomes larger as the strength of the interaction keeps 
increasing. 



*Elcctronic address: jianpinwu@mail.bnu.cdu.cn 
^Electronic address: zcnghbi@gmail.com 



1 



I. INTRODUCTION 



As a new method to study strongly coupled quantum field theory, AdS/CFT correspon- 
dence [l-3] has been used to study the strongly coupled phenomena in quantum many-body 
systems. One of the most successful progress is that by coupling a free Fermion field to the 
Reissner — Nordstrom (RN) black hole, a class of non-Fermi liquid at zero temperature 
was found holographically 4j-[6j • The properties of this non- Fermi liquid in the presence of 
a magnetic field have been studied in 7H9(. The extension to finite temperature was investi- 



gated in 



101 ] . In ll|, |l2j, free fermions in AdS BTZ black hole and Gauss-Bonnet black hole 



have also been studiec 
been explored in Refs. 
on this subjects. 



respectively. In addition, the fermions on Lifshitz Background have 



13 



141 ] . We can refer to the lecture [15| for an excellent latest review 



In 



16 



171 ]. the authors studied the holographic fermion system with RN black hole 



background when a dipole interaction is added, it is found that the strength of the interaction 
p is similar to U/t in the fermions Hubbard model 1 . When p is large enough, a gap opens just 
like in a Mott insulator, in which a large repulsive interaction U will change the structure of 
the bands. Then the fermion energy is in the middle of the gap, which makes an insulator 2 . 

However, the RN black hole background has nonzero ground state entropy density, which 
seems to be inconsistent with our intuition that a system of degenerate fermions has a unique 
ground state. Therefore, a systematic exploration of the system that has zero extremal 



entropy will be important and valuable. Such models have been proposed in Refs. 



21 



231 ] . Furthermore, in Ref. 24j, the author investigated the free fermionic response in the 



background proposed by Gubser and Rocha 



211 ] . They find that the dispersion relation is 



linear, just like a Fermi liquid. It is very different from that found in RN black hole 



Among these models of zero ground state entropy density, a very important and valuable 



model is that proposed by Goldstein et al. 



221 ] . Here the near horizon geometry is Lifshitz- 



like. Therefore, in this paper, we will study the characteristics of the fermionic response 



1 InRefs.0,0, the authors construct the Pauli couplings in a bottom-up setup. The inspiration for using 
Pauli couplings in holographic AdS / CMT came from the previous work 18|, |l9( , in which they derived the 
type of fermionic interactions that appear in the consistent siiriergravity truncations. 

2 Another important paper on the effect dipole coupling is Ref. 2fJ , which mainly focus on the smaller value 
of dipole coupling. They find that the non-Fermi liquids behavior is robust under this deformation. 
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with dipole interaction in this background 3 . The interesting results we find is that, similar 
to the fermions in the RN black hole background, the dipole interaction will open a gap for 
the strongly coupled fermions system when the value of the interaction p is large enough. 

The organization of this paper is as follows. In section II, we first review the charged 
dilaton black branes with zero entropy at zero temperature. In section III, we derive the 
equations of motion for the fermion field with a dipole action. The main results of the 
emergence of gap is in section IV. Finally, the discussion and conclusions are presented in 
section V. 



II. AN EXTREMAL CHARGED DILATONIC BLACK BRANE SOLUTIONS 



A. Einstein-Maxwell-dilaton model 



For completeness let us also review the dilaton gravity system 22] . Our discussion follows 



22i | closely and some details are omitted. The system we consider is as follows 4 : 

S = ^J d 4 x^[R - 2{d<P) 2 - e 2 ^F ab F ab + A], (1) 

where R is Ricci scalar, is the dilaton field, F a b = d a Ab — d b A a is the field strength, and L 
is the AdS radius. In this paper we will focus on an electrically charged black brane solution. 
The metric and gauge fields can be written as, 

dv^ 

ds z = - a \r)dt 2 + + b 2 (r)[(dx) 2 + {dy)\ (2) 

a (r) 

e 2a<t>(r) F = J^dt A dT. (3) 

tr (r) 

Such a gauge field satisfies the gauge field equation of motion automatically. With the above 
ansatz, the remaining equations of motion can be expressed as follow: 

(a 2 &V)' = -ae" 2 ^, (4) 



3 The free fermionic response in the present background and the fermionic response with dipole coupling 
in the background propose d hv Clii hser and Rocha, will be addressed in two companion papers. 

4 We can also refer to Refs.[25|, [26j, in which the holographic models of charged dilatonic black branes 
(without charged fluid) at finite temperature and several classes of coupling functions are discussed in 
some details. 
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a 2 b' 2 + Ua 2 )'(b 2 )> = <pW - e~^%- + %. (5) 

lib 2 

W = ~§r- (6) 
h" 

T = ~ (7) 

where the prime denote the differentiation with respect to the coordinate r. In the subse- 
quent subsections, we will construct such solutions that the near-horizon geometry has a 
Lifshitz-like symmetry but at infinity the geometry asymptotes to AdS. For convenience, we 
shall set L — 1 and k — 1 in the following. 



B. Scaling solution near the horizon 

In order to construct the scaling solution near the horizon r^, we consider the following 
ansatz: 

a = drl b = C 2 r?,</> = -K In r* + C 3 , (8) 

where C±, C 2 , C 3 , 7, f3, and K are all constants and C\ and C 2 are set to be positive in 
the following, without loss of generality. Note that for convenience, one has introduced the 
variable r* = r — in the above ansatz. With this ansatz and the equations of motion, we 
have 



6 o2 aoca = (2/3 + lj^gg 

l + (f) 2 '" l + d) 2 '" 1 (/3 + l)(2/3 + l)' V « 

(9) 



Such a solution has a Lifshitz-like symmetry in the metric, with a dynamical critical 
exponent z = -s, although such a symmetry is broken by the logarithmic dependence of the 
dilaton on r*. Also we note that in the above solution, the metric component gu has a double 
zero at the horizon where g xx also vanishes and therefore corresponds to an extremal brane 
with vanishing ground state entropy density. Finally, we also point out that for a = 0, such 
a solution corresponds to an AdS% x R 2 geometry, which is also the near horizon geometry 
of the extremal RN black hole. The fermionic response in RN black hole has been explore 
carefully in Refs.j^j, 6]. Therefore we shall assume a > and focus on the Lifshitz like near 
horizon geometry in the following. 
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C. The solution with asymptotes to AdS at infinity 



To obtain a solution which asymptotes to AdS at infinity, we can add a perturbation to the 
scaling solution which is irrelevant in the IR scaling region but relevant in the UV region. 
In order to set the constant C 2 to unity and C 3 to vanishing, we carry out a coordinate 
transformation under which, 

r = \r,t = ^,x i = ?L (10) 
A A 

c 3 



with A = e k and A = \[C 2 }? ■ Under such rescaling, we have 

a = Cm,b = r^,<j) = -Klnn, (11) 
and in terms of a, the charge can be expressed as 

< 12 > 

Subsequently, we add the perturbation to the Eqs. (jSJ). The resulting functions satisfied 
the equations of motion at the leading order are 

a = dr*(l + diO, b = + d 2 <), = -iHnr* + d 3 <, (13) 

where Ci, /3, and A' keep unchanged. The perturbation is characterized by the exponent 
v and the three constants d\, d 2 and d 3 . Also, c?2 and c?3 are determined by d\ and they 
have the following relations: d 3 = 2/3 ^-~ 1 (j 2 , c?i = [ ^ 2 p+^v)(2p+i+v) ~ In addition, the 
perturbation should die out at small r*. Therefore, we must require v > 0, which gives rise 



to a unique allowed value for this exponent: v = |[— 2/3 - 1 + ^(2/3 + l)(10/3 + 9)]. Thus 
there are two parameters a and d\ for this perturbated solution. For simplicity, we will set 
a — 1. 

The initial data for numerical integration is taken from the above perturbated solution 
near the horizon. As depicted in Fig. (TJ for the perturbation with d\ = —0.514219, b'(r) — > 1 
as r — > oo, that is to say, the solution near the boundary takes the standard asymptotics of 
AdS for the perturbation with d\ = —0.514219 5 . In the next section, we shall investigate 
the holographic fermions with dipole term in such a background. 



5 For any other negative value of d\ , b does not take the standard form of asymptotic AdS and we need a 
coordinate transformation to achieve our goal. When d\ > 0, the numerical solution is singular. 
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FIG. 1: Left Plot: Background solution in the case of a = 1 and d\ = —0.514219. Blue line denote 
a(r), red line denote b(r) and green line denote 4>(r). Right plot: The plots for a'(r) (blue line) 
and b'(r) (red line). It indicates that a'{r) and b'(r) approach 1 when r approaches oo. 

III. HOLOGRAPHIC FERMION WITH BULK DIPOLE COUPLING 



In order to know the effect of magnetic dipole coupling on the structure of the spectral 
functions of the dual fermionic operator, we consider the following bulk fermion action 6 



s D = % J d 4 xV^c (ra a 



m 



ipft) C, 



(14) 



where T a is related to the usual flat space gamma matrix by a factor of the vielbein, T a = 
(e^) a r At , V a = d a + ^{uj^aY^ — iqA a is the covariant derivative with {u^a the spin 
connection 1-forms and ft = ^T fJiU (e fl ) a (e u ) b F a i J . Before taking the particular background 
(Eqs. (j2J) and (j3J)), we will derive the flow equation as Refs.[4 
static background, i.e., 



12 



24] in a more general 



-g tt (r)dt 2 + g rr (r)dr 2 + g xx {r)dx 2 + g yy (r)dy 2 , A a = A t (r)(dt) a . 



Firstly, the Dirac equation can be derived from the action Sd 

T a V a C - mC - ipft( = 0. 



(15) 



(16) 



2, 2 it 

Since the gauge field A t can be expressed as A t = a qJ (It can easily be obtained by Eqs. © and (gj)), 
the bulk fermion action ()14|) . only through the combination of (qQ,pQ), depends on (q,p). In this paper, 
for convenience, we only focus on the case of q = 4 and p > 0. For more general parameter space 
we will discusses them in the near future. When we have set a = 1, Q = ±y/2. Without loss of generality, 
we shall take Q = \/2 in the following. In this case, the chemical potential can be obtained as /i ~ —0.843 
numerically. If we take Q = — \/2, the same results will be recovered by tuning the parameter (q,p). 
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Making a transformation £ = (—gg rr ) 3 J 7 to remove the spin connection and expanding 
J 7 as J 7 = p Q-^t+ikiX 1 j n p our j er space, the Dirac equation ( fl6l) turns out to be 



(v/7rr r d r - m - —^g rr g tt T rt d r A t )F - i(co + qA^^^F + ik^g^Y x F = 0, (17) 

where due to rotational symmetry in the spatial directions, we set k x = k and ki = 0, i ^ x 
without losing generality. Notice that Eq. ffT7|) only depends on three Gamma matrices 
T r ,T t ,T x . So it is convenient to split the spinors F into F = (Fi,i^) T and choose the 
following basis for our gamma matrices as in 6| : 



-a 3 
-a 3 



to 1 
ia l 



-a 2 
a 2 



(18) 



So, we have a new version of the Dirac equation as 



fg^dr + 777.(7 



F 2 



g tt (oj + gA)ia 2 =F kyfg^a 1 - p\/g tt g rr d r A t a 1 



Fx 



.(19) 



Furthermore, according to eigenvalues of T r , we make such a decomposition F± = |(1 ± 
T r )F. Then 



Si 
£ 2 



^2 



(20) 



Under such decomposition, the Dirac equation ffT9l can be rewritten as 



\ Al ) =±(u + qA t )y/<p \\- (VF * + p\fg tt g rr d r A t ) [ Bl 
Bt UJ Ui 



(21) 



fg^d r ± m) | = ±(u + g^) j + (k^ - p^g^d r A t ) ^ | . (22) 

Introducing the ratios £ a = = 1,2, one can package the Dirac equation (I21j) and 



22]) into the evolution equation of £ a , 



/^<9 r + 2m)£ Q = [v. + (-l) a AV<F] + [i 



-i) a kVg^] e 



2 

O 5 



(23) 



where t>± = y/g^{uj + gA 4 ) ±pyf g tt g rr d r A t . The above flow equation will be more convenient 
to impose the boundary conditions at the horizon and read off the boundary Green functions. 



Now, we take the particular background (Eqs. (|2J) and (151)). Because near the boundary, 
the geometry is an AdS^, the solution of the Dirac equation (fT9]) can be expressed as 

F a r ^°a a r m 1° J +b a r~ m f* J , a = l,2. (24) 



are related by 6 C 











If 6« 1 




and a a 













functions G is given by G = —iS^° 29|. Therefore 











= Sa a | 











G(w, jfe) = lim r 2m I ^ ° I , (25) 



£ 



2 



Near the horizon, we take the scaling solution near the horizon (jSj). We find that the 
requirement that the solutions of Eqs. (1211) and (122!) near the horizon be in-falling implies 

&/ = H /or w^O. (26) 



IV. EMERGENCE OF THE GAP 



In this section, to study the spectral function, we will numerically solve the flow equation 
(123]) with the initial condition ( 126]) . Up to normalization, the spectral function is given by 
A(oj, k) = Im[TrGji(u, k)]. However, due to the relation G u (u, k) = G 2 2(w, —k), Gu(u>, k) 
can be recovered from G^w, k). Therefore, in the following 3D plots and density plots, to 
make numerical calculation easier, we will only show G 22 (a;, k). 



A. The Fermi momentum kp 



Before discussing the spectral function for the larger p, we would like to give a simple 
discussion on the Fermi momentum kp for p = 0. From the plots above in FIGJ2] one can 
see that for p = 0, a sharp quasiparticle-like peak occurs near k ~ 0.75 and cu = 0, indicating 
a Fermi surface. When q is increased, the Fermi momentum kp increases (Table I). As q is 
further increased, several Fermi surfaces emerge. For example, for q — 6, the another Fermi 
surface begin to occur at k ~ 0.168. In fact, the case that with the increase of q, several 
Fermi surfaces produce also occurs in RN-AdS background. It seems that the gravity duals 
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FIG. 2: The 3d and density plots of I mG 22^1 k) for q = 4 and m = (the plots above for p = 
and the plots below for p = 2.). In the plots above, a quasiparticle-like peak occurs at k ~ 0.75 
and w = 0. However, in the plots below, a gap emerges around uj = 0. 

with larger charge g posses more branches of Fermi surfaces. Conversely, with the decrease 
of q, the Fermi sea gradually disappears. 
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Q 


3 


3.5 


4 


4.5 


5 


5.5 


6 




0.2481 


0.4874 


0.7565 


1.0457 


1.3494 


1.6643 


1.9878 



TABLE I: The Fermi momentum kp for different charge q. 



In addition, we would also like to point out that for fixed q — — 1 and a = 10 -6 , hp ~ 
0.917, which recover the results of Ref.j^] 7 . For more discussions on the Fermi momentum 
kp and the low energy behaviors in this case, we will explore it in the another companion 
paper. 



B. Emergence of the gap 



When we turn on p to the value of 2, a gap emerges and there are two bands located 
at the positive frequency and negative frequency regions, respectively (the plots below in 
FIGJ2]). Evidently, the strength of the lower band is bigger than the upper band. We also 
show the 3D plot and density plot of the Green's function lmG<i2 for p = 2.5 and p = 4.5 
in FIGJH We can see that with the increase of p, the width of the gap becomes larger and 
the band in the negative frequency region switch gradually to the positive frequency region. 
As p increases further, the lower band disperses and the upper band becomes stronger and 
sharper. 

In order to explore further the characteristics of the spectral function, We show the 
spectral function A(u, k) for p = 0, 2, 2.5 and 4.5 for sample values of k. From the left plot 
above in FIGJH we find that some peaks distribute at both positive and negative frequency 
regions. When the peak approaches u = 0, its height goes to infinity, and its width goes to 
zero, indicating that there is a Fermi surface at k = kp. When we amplify p, the strength 
of the quasiparticle-like peak at u = degrades, and vanishes at a critical value p cr u- For 
p > p cr u, the dipole interaction will open a gap as the case in RN black hole 16|. Differ from 
the case of p < p cr it, the height of the peaks degrade and the width becomes larger when 
u> = is approached, and the peak vanishes around uo = 0. We also note that the spectral 
density mainly appears at negative frequency regions for p = 2 (the right plot above in 



7 Since we take the negative chemical potential here, q = — 1 corresponds to q = 1 in Ref. [J|. In addition, for 
a = 0, the numerical solution is singular. Therefore, in the numerical calculation, we take a = 1CP 6 ~ 
instead of a = 0. 
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FIG. 3: The 3d and density plots of IraQiii^^ k) for q = 4 and m = (the plots above for p = 2.5 
and the plots below for p = 4.5.). 

FIGS]). As p increases, the spectral density begins to distribute the positive frequency and 
negative frequency regions, respectively (the left plot below in FIGHJ). When p increases 
further, the spectral density switches gradually to the positive frequency regions (the right 
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FIG. 4: The plots of spectral function A(uj, k) as a function of oj for sample values of A; € [0.6, 3.8] 
for p = (left plot above), p = 2 (right plot above), p = 2.5 (left plot below) and p = 4.5 (right plot 
below). For p = 0, the sharp quasiparticle-like peak occurs at u) = when we dial k = kp ~ 0.75. 
However, for p = 2, 2.5 and 4.5, the gap around u = persists for all A;. 

plot below in FIGJ4J) 8 . 

The another important quantity of interest for us is the density of states A(u), which is 
defined as the integral of A(u), k) over k. It is the total spectral weight. By careful numerical 
computations, we find that the onset of the gap is near p = 1.5, which is different from the 
onset value of p = 4 in the case of RN background. Consistent with the above observations 
(FIGfSl FIG J21 and FIGHJ), the total spectral weight is also mainly distribute at negative 
frequency regions for small p (p > p cr it, left plot in FIGJ6]), and switches to the positive 
frequency regions as the increase of p. 

8 It seems that for the case of the RN black hole, the spectral density always appears at negative frequency 
regions (p > 0,FIG[5]). 
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A(w,k) 




A(u,k) 




FIG. 5: For the background of RN black hole, the plots of spectral function A(lj, k) as a function 
of ui for sample values of k £ [1, 4.5] for p = 6 (left plot), p = 8 (right plot). 
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FIG. 6: Left plot: The density of states A(u>) for p = 1 (green), p = 1.5 (blue) and p = 2 (red); 
Right plot: The density of states A[u) for p = 2.5 (blue), p = 3 (green), p = 4 (red) and p = 4.5 
(black). The onset of the gap is at p ~ 1.5. 

Finally, we show the relation between the width of the gap A and p. Evidently, the gap 
becomes wider as p increases. 

V. CONCLUSIONS AND DISCUSSION 

In this paper, we have explored the fermionic response in the presence of a bulk dipole 
interaction term with strength p in the dilatonic black brane with a Lifshitz like IR geometry 
and AdS± boundary. When p = 0, a sharp quasiparticle-like peak occurs near k = 0.75 and 
to = 0, indicating a Fermi surface. As p increases, the peak become wide and its strength 
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FIG. 7: The width A of the gap as a function of p. 



degrades gradually, and the spectral weight is transferred between bands. Furthermore, 
when p goes beyond a critical dipole interaction p cr it, the Fermi sea disappears and a gap 
emerges as the case in RN black hole 4j. By studying the another important quantity: the 
density of states A(u), we clearly find that the onset of the gap is near p = 1.5 and the 
width of the gap becomes bigger as p increases. 

So far, we have found a (Mott) gap from holographic fermions in the background of the 
dilatonic black brane with a Lifshitz like IR geometry when a dipole interaction term are 
added as that found in RN black hole. However, for more general dilatonic black brane, 



which a 
in Refs. 



so have vanishing ground entropy density, for example, a class of systems proposed 



21 



32 



34|, is the emergence of gap robust? The fermionic response without dipole 



interaction in this class backgrounds has been investigated in Refs. 24j, |35|. Therefore, it 
is valuable to test the robustness of the emergence of the gap in the backgrounds of more 
general dilatonic black branes. In addition, it is also interesting to extend our investigations 
in this paper to the case of the non-relativistic fermionic fixed point. We will address them 
in the near future. 
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